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HIGHLIGHTS 


• We analyze heat transfer enhancement in a liquid-piston compressor for CAES. 

• Hydrothermal characterizations of the inserted heat exchanger media are obtained. 

• CFD simulations are done combining the VOF and porous media modeling methods. 

• The heat exchangers can effectively suppress temperature rise and secondary flows. 

• Results indicate optimization on the shape of the heat exchanger could be done. 
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In this paper, the topic of Compressed Air Energy Storage (CAES) is discussed and a program in which it is 
being applied to a wind turbine system for leveling power supplied to the grid is described. Noted is the 
importance of heat transfer in the design of the compressor and its effect on performance. Presented is a 
design for minimizing the temperature rise in the compressor during compression. The design requires 
modeling regenerative heat transfer from the compressed air to solid material inserted in the 
compression space. Modeling requires characterizing pressure drop through the porous insert, interfacial 
heat transfer between solid and fluid in the matrix, and thermal dispersion within the porous regions. 
Computation and experimentation are applied for developing correlations for such terms. Two types of 
porous media are applied: interrupted plates and open-cell metal foams. Cases with foam inserts are 
computed and the results are discussed. Discovered in the results are some complex secondary flow 
features in spaces above the porous inserts. 

© 2014 Elsevier Ltd. All rights reserved. 


1. Introduction 

This paper presents thermal analyses on a liquid piston driven 
compressor used for Compressed Air Energy Storage (CAES). The 
CAES system stores energy as high-pressure air, to retrieve it later 
in a liquid piston expander. Compression leads to a tendency for 
temperature rise in a compressible gas. Absorbing heat from air 
during compression to reduce its temperature rise is important for 
improving compression efficiency. As the air temperature rises, 
part of the input work is being converted into internal energy rise 
that is wasted during the storage period as the compressed air cools 
toward the ambient temperature. This paper discusses techniques 
for minimizing that temperature rise. 
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Our desire is that the CAES system is integrated into a wind 
power field. The benefit of integrating energy storage into a wind 
power array is that fluctuations in wind power input can be 
smoothed over time and electric power generation equipment can 
be sized commensurate with a supply power that is nearer the 
average wind power of the day. The alternative is sizing for peak 
power or “throttling back” the wind turbine when the wind is 
strong. 

Use of storage techniques for wind power is discussed in Ref. 1 . 
Various scenarios with and without storage were considered. In 
one, wind generation capacity in 2050 was 302 GW without storage 
and 351 GW with storage. Another study, which is modeled using 
CAES and wind energy to supply 20% of the total power generation 
in Texas concluded that the addition of CAES to wind energy re¬ 
duces the impact of wind fluctuations on the electricity grid [2]. 
From a report written by the US Offshore Wind Collaborative [3], 
we read that the total amount of U.S. offshore wind energy capacity 
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Nomenclature 

a v surface area per unit volume of porous medium (1/m) 

A cross-sectional area (m 2 ) 

b half plate distance (m) 

/ coefficient for the Forchheimer term (1/m) 

c p constant-pressure specific heat (J/kg I<) 

c s specific heat for solid (J/kg I<) 

c v constant-volume specific heat (J/kg K) 

Dh hydraulic diameter (m) 

D characteristic diameter (m) 

d characteristic length based on filament dia. (m) 

df filament diameter (m) 

d m mean pore diameter (m) 

E s storage energy (J) 

g gravitational acceleration (m/s 2 ) 

/i s f surface heat transfer coefficient (W/m 2 K) 

hy volumetric heat transfer coefficient (W/m 3 I<) 

I< permeability (m 2 ) 

k thermal conductivity (W/m K) 

/< dis dispersion conductivity (W/m K) 

L chamber length (m) 

L\ length of the upper region without insert (m) 

L ins length of the insert region (m) 

l plate length (m) 

m total mass of air (kg) 

n polytropic exponent 

Nu Nusselt number 

P bulk pressure (Pa) 

Ptot supplied hydraulic pressure (Pa) 

Pe Peclet number 

q area averaged wall heat flux 

p local pressure (Pa) 

Pf final air pressure (Pa) 

Pr Prandtl number 

R radius of chamber (m) 

Pei Reynolds number based on 1 

r radial coordinate 

ideal gas constant (J/kg K) 


Re Reynolds number based on characteristic length d 

S m momentum source term (Pa/m) 

T local air temperature (K) 

T 0 initial temperature; wall temperature (K) 

T a average air temperature in the chamber (K) 

T s local solid temperature (K) 

Psoiid average temperature of solid in the chamber (K) 

T f final averaged air temperature (K) 

t plate thickness (m) 

/ time (s) 

U [n liquid piston velocity (m/s) 

V instantaneous volume of chamber (m 3 ) 

W in work input (J) 

x axial coordinate 

Greek symbols 
a volume fraction 

£ porosity 

7] compression efficiency 

/j, dynamic viscosity (Ns/m 2 ) 

p density (kg/m 3 ) 

£ final pressure compression ratio 

Subscripts 

0 initial value of variable 

1 air phase 

2 water phase 

D based on filament diameter 

Da Darcian 

f fluid phase 

REV averaged on the REV 

s solid 

XX local streamwise direction 

YY local cross-stream direction 

Superscripts 

* dimensionless variable 


is almost equal to the current total installed capacity. Noting that 
most electricity-demanding regions are on the coasts, the report 
highlights the significant potential for offshore wind generation. 
Traditional CAES systems store the compressed air in underground 
caverns and often marry this with combustion plants by using 
compressed air as oxidant for burning in combustion turbines [4]. 
More recent advancements include using multi-stage compressors 
with inter-cooling between stages and multi-stage expansion tur¬ 
bines, with reheat [5 . A small-scale wind turbine system with CAES 
was constructed and tested [6 . Although the study showed bene¬ 
ficial prospects of using combined CAES with the wind turbine 
system for home applications, the efficiency of the system tested 
was very low, mainly due to a lack of cooling during compression. 
Herein, near-isothermal compression and expansion are by using 
regenerative heat exchange with elements within the compressor/ 
expander volumes. 

In the present system, we employ the “open accumulator” 
concept proposed in Ref. [7 . In this system, gas is exhausted to the 
atmosphere during expansion, reversing during compression. 
Compression is with pumped liquid and air and expansion is with a 
liquid motor and an air expander. It is operated at a constant 


accumulator pressure; the maximum design pressure. Thus, with 
the open accumulator, the energy storage density per unit volume 
of air is always high since the low-power, low-accumulator-pres¬ 
sure situation in the closed accumulator is avoided. An added 
benefit is a dramatic reduction in accumulator pressure oscillation 
cycles and improved fatigue performance. 

Though not necessary for the open accumulator CAES to be 
successful, our design is proposed to be used in a fluid power wind 
turbine system. In such as system, we use hydraulic equipment for 
power conversion and transmission, as shown in Fig. 1 . It employs 
hydraulic circuits, a hydraulic drive pump, hydraulic pumps and 
motors, open accumulators, and liquid-piston air compressors [8]. 
Advantages include having a hydraulic pump in the nacelle at the 
top of the tower rather than a heavy transmission and generator 
assembly and having much of the equipment, including the 
generator, at ground (or sea) level, allowing easier access for as¬ 
sembly and maintenance. When electricity demand is small, the 
system operates in storage mode in which the excess shaft power 
from the wind turbine is transferred to hydraulic power in the drive 
pump and the hydraulic power is transmitted through the hy¬ 
draulic circuit to the hydraulic transformers, which are a series of 
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Fig. 1. Schematic of CAES system for offshore wind energy storage and generation [8]. 


specialized hydraulic pumps and valves to pump the liquid into the 
liquid-piston chamber to compress air. During this mode, com¬ 
pressed air is stored in the open accumulator. When the electricity 
demand is larger than the direct supply from the wind turbine, the 
unit switches to the generation mode in which air is withdrawn 
from the open accumulator and expanded in the liquid piston 
chamber to derive shaft power for the generator. 

Key to success of the system is efficient compression and 
expansion of air in the compressor/expander. As noted, this must be 
done under near-isothermal conditions. Thus, the present paper 
will concentrate on the thermal design of the compressor/expander 
and the analysis on which it is based. This paper focuses on 
compression and expansion with a liquid piston. The term “liquid 
piston” is applied because the gas compression is done with a liq¬ 
uid-gas interface rather than a solid piston. It is found in Ref. [9] 
that liquid-piston compressors have an advantage over solid pis¬ 
tons in terms of efficiency and power consumption. Further, for the 
purpose of cooling, insertion of solid material into the chamber is 
possible with liquid-piston compressors as liquid can flow through 
the open portions of the matrix. The use of porous media and 
regenerative heat exchangers to absorb thermal energy has been 
shown to be successful in a number of applications, such as a micro- 
fabricated, segmented, involute-foil regenerator used in a Stirling 
engine 10], a honeycomb regenerator used in a high-temperature- 
air-combustion system 11 , a regenerative heat exchanger used in 
a reverse-flow reactor 12], and a magnetic regenerative system 
used in a refrigeration cycle 13]. In the present study, it will be 
shown that having a liquid piston design for the compressor 
inserted with porous media is a major step in having an efficient 
compressor/expander. The remainder of the paper will be dedi¬ 
cated to modeling thermal storage to the porous material in the 
compressor volume in preparation for design of the compressor. 

2. Analysis on heat transfer and compression efficiency 

The effect of heat transfer on compression efficiency will be 
shown quantitatively through a simple thermodynamic analysis. 


The compression efficiency for CAES systems is defined as the ratio 
of storage energy to work input [14,15 . The storage energy is 
defined as the amount of work extraction from compressed air as it 
is isothermally expanded to atmospheric pressure. 

E s = mm 0 (in P* - 1 + T) (1) 

Compression of air is completed in two steps: compress the air 
from atmospheric pressure to a high pressure, typically resulting in 
temperature increase, and allow the air to cool while compressing 
at the storage pressure so that the work potential (storage energy) 
is maintained as volume decreases. Work of the latter stage is 
identified as “cooling work.” The total work input is the sum of the 
work done in the two steps. It is given by. 

v f 

Win = - J (P - Po)dV + (Pf - P 0 ) (Vf - (2) 

Vo 

The compression efficiency is given by the ratio of the storage 
energy to the total work input. 


V 


mm 0 (in P* - 1 + A 


W-, 


in 


(P-P 0 )dV+(P 0 K-P 0 )(Vf- Vo 


( 3 ) 


V 0 


K 


It is shown that the compression process can be characterized as 
a polytropic process [15 , 


PV n = P 0 V£ (4) 

where the polytropic exponent, n, shows the effect of heat transfer; 
n = 1 represents isothermal compression and n = 1.4 for air rep¬ 
resents adiabatic compression, with no heat transfer. Substituting 
Eq. (4) into Eq. (3), and writing in terms of dimensionless pressure 
and volume based on their initial values, 


lnP*-l+i 

V = - (5) 

p,( ' 1 '^ / 1 "- 1 + P’-Vn - 1 + (P* - 1) (p*-l/n - 1) 

Eq. (5) shows that the compression efficiency is dependent on 
the pressure compression ratio, P* and heat transfer, given by n. 
Fig. 2 shows that for the same pressure compression ratio, 
decreasing n, which means enhancing heat transfer from the 
compressed air, significantly improves efficiency. 



Fig. 2. Effects of heat transfer on compression efficiency with different pressure 
compression ratios. 
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3. Heat transfer to the porous inserts 

3.1. Introduction to porous inserts 

Porous media can be inserted into the liquid-piston compres¬ 
sion chamber so that the solid can absorb thermal energy from the 
air as it heats during compression. The present study will introduce 
two kinds of porous inserts, an open-cell metal foam and an 
interrupted-plate insert (Fig. 3). The metal foam has very thin 
filament features that enhance mixing of the flow. The interrupted- 
plate insert consists of an array of plates oriented in a staggered 
fashion so that a new thermal boundary layer develops on each 
successive plate in the streamwise direction. Both have large heat 
transfer surface area to volume ratios. The open-cell metal foam is a 
common type of regenerative heat exchanger. The interrupted- 
plate insert provides opportunities for us to design its shape to 
result in improvements; the form of the matrix is chosen, but then 
different sizes of the shape parameters are studied to adjust the 
shape. 

The flow resistance and resulting pressure drop, the heat 
transfer coefficient between the solid material and surrounding air, 
and thermal dispersion will be discussed for the two types of 
porous inserts discussed in this paper. Then, using these models as 
source terms in the transport equations, simulations on the full 
compressor with a metal foam insert will be done. 

3.2. Interrupted plate 

Because details of the shape of the interrupted-plate insert are 
known, CFD simulations on a Representative Elementary Volume 



(a) Metal foam inserts 



(b) Interrupted-plate insert 



(c) A unit cell of the interrupted-plate insert 


(REV) of the interrupted-plate can be done. These CFD analyses are 
used to obtain the pressure drop and interfacial heat transfer re¬ 
lationships for flow with various velocities through the matrix. The 
REV is the minimum representative, repeating geometric feature of 
the porous medium. Its geometry for the interrupted plate insert is 
shown in Fig. 3(c). Three geometric dimensions are: the plate 
length, /, the plate distance, 2b, and the plate thickness, t. 


3.2.1. Shape analysis 

A study was conducted to quantify the effects of the shape pa¬ 
rameters on heat transfer and compression efficiency 16]. Twenty- 
seven different shapes were analyzed. For each shape, two flow 
conditions were studied: Re\ = 1 and Re\ = 400, where Re\ is defined 
as. 


P^REV,f 

MREV,f 



The flow is presumed to be in the same direction as the x axis 
(Fig. 3(c)). Periodic momentum and thermal boundary conditions 
are used on the entry and exit surfaces of the REV model such that 
the flow being simulated is representative of a fully-developed flow 
through the central region of the porous matrix. The solid surfaces 
are maintained at a uniform and constant temperature in the REV 
simulations. 

From these CFD simulation results, a heat transfer correlation, 
generalized in terms of hydraulic diameter, is developed. The 
Reynolds number and Nusselt number based on hydraulic diameter 
are: 



pPh u REVJ 

MREV,f 



Nu D h = (8) 

where, 

h = =- \ - (9) 

h ~ 1 bulk,f 


Dh = 4b (10) 

The heat transfer correlation is given by Ref. [16]: 

Nu Dh = 9.700 + 0.0876Re^ 92 Pr3 (11) 

It is plotted with the data from the various CFD runs in Fig. 4(a). 
As the shape parameters in the CFD runs can be varied, this cor¬ 
relation captures the effects of different shapes on heat transfer. 

The pressure drop for flow passing through the interrupted- 
plate exchanger can be modeled using a viscous and an inertial 
term. 


dp p F 2 

0^ = -^UREV - PUrev 


( 12 ) 


Using data from the study given in Refs. [16 , the permeability 
and Forchheimer coefficient for a given plate geometry is correlated 
against the plate length and plate thickness by a least square fit: 



0.1151 lg( —) — 5.184^--173.9^— + 1.031 

ylm ) lm lm 

x10 -6 m 2 


Fig. 3. Porous inserts used in the liquid piston compressor. 


(13) 


Please cite this article in press as: C. Zhang, et al., Thermal analysis of a compressor for application to Compressed Air Energy Storage, Applied 
Thermal Engineering (2014), http://dx.doi.Org/10.1016/j.applthermaleng.2014.08.014 























ARTICLE IN PRESS 


C. Zhang et al / Applied Thermal Engineering xxx (2014) 1—10 


60 

50 

40 

130 
20 
10 
0 


Nu Dh = 9.700 + 0.0876/?<? nh ° 792 PrI > 


Dh 


• -Fit 


♦/ 

/ 

/ 


♦/ 

/ 




a Different 
ReJ 

♦ Re_l=l 

♦ Re 1=400 


0.1 1 10 100 1000 

R e Dh 

(a) Dimensionless numbers based on hydraulic diameter, D h , 
computed from CFD runs on the REV with various values of the 
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(c) Dimensionless heat transfer correlation for shape: l = 7.5 mm, 

t = 0.8mm, 2b = 2.5mm 
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(b) Area-averaged wall heat transfer coefficient for shape: / = 7.5mm, 

t = 0.8mm and 2b = 2.5mm 
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Fig. 4. Various heat transfer simulation studies on the REV of the interrupted-plate exchanger. 


1.057 


F 2 . 5 = -0.009979 lg ( +j + 62.64(7 ] - 0.04957 


lm 


lm 


(14) 


I<5 


0.3449 lg( 7-) + 1.285—-— 394.0 7— + 3.195 
\lm ) lm lm 

xlO -6 m 2 


(15) 


F 5 = 0.003209 lg(7-) - 1.1347- +403.0(7- 

\lmj lm \lm 

+ 0.03676 


1.493 


(16) 


IC 


10 


— 0.03438 lg( zr—~] + 155.3 1 


lm 


lm 


984.2+ + 3.835 
lm 


x 10 6 m 2 


(17) 


F 10 = 0.009688 lg( -—) - 1.319++ 9.397+ + 0.06492 

ylm ) lm lm 

(18) 

where subscripts 2.5, 5 and 10 represent plate lengths of 2.5 mm, 
5 mm and 10 mm, respectively. 

Using the interfacial heat transfer correlation, the compression 
process is simulated based on a zero-dimensional, thermodynamic 
model, from which the efficiency (Eq. (3)) is evaluated for each 
shape [16 . It is found that, in general, decreasing the separation 


distance between plates in the interrupted plate insert leads to 
improved efficiency. For a given plate separation distance, 
decreasing the plate length improves efficiency, as shown from a 
simple boundary layer development analysis. The heat transfer 
correlation and the pressure drop equations developed from this 
study can be used also in optimization of the exchanger shape. 


3.2.2. Heat transfer and pressure drop under high Reynolds number 

Analyses in this section are focused on heat transfer under 
relatively high Reynolds number flow. Simulations are done at 
eleven different Reynolds numbers, Re D h, ranging from 0.67 to 
5333, for an interrupted-plate insert with the following shape pa¬ 
rameters: l = 7.5 mm, t = 0.8 mm, 2b = 2.5 mm. Three different runs 
are computed at Re D h = 5333 while varying the density of the air 
flow and the mass flow rate. Each is a single CFD run. The com¬ 
mercial CFD Software ANSYS FLUENT is used and the k-e model is 
used for turbulence closure. The REV models have periodic mo¬ 
mentum and thermal boundary conditions. Details of implement¬ 
ing these conditions are discussed in Ref. [16 . The number of 
computational cells varies with Reynolds number; specifically: 
169,632 for Re Dh < 150,588,816 for 150 < Re Dh < 2000,1,357,056 for 
2000 < Re oh < 5000, and 2,024,352 for Re D h > 5000. These grid 
cells have been verified by grid-independence studies. From the 
CFD simulations, heat transfer coefficients and pressure drop terms 
are obtained. 

The area-averaged surface heat transfer coefficient depends on 
both flow velocity and density, as shown in Fig. 4(b). The Reynolds 
numbers and Nusselt numbers can very well capture heat transfer 
under different velocity and pressure conditions, as shown in 
Fig. 4(c). The permeability and Forchheimer coefficient for this 
shape are also computed: 
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I< = 2.813 x 10“ 7 m 2 , F = 0.0359 

Another set of simulations are done at four Reynolds numbers 
for an interrupted-plate insert with shape parameters: l = 7.5 mm, 
t = 0.55 mm, 2b = 2.75 mm. For Re oh = 1.3 and 286.9, a grid of 
707,840 cells is used and the laminar flow is solved. For 
Re Dh = 13,415.6 and 66,820.6, the grid cells are, respectively, 
1,469,139 and 2,584,833, and the Transition SST model (see 
Ref. 17]) is used for turbulence modeling. The Nusselt numbers and 
Reynolds numbers of these simulations are compared with results 
from all previous simulations on the basis of hydraulic diameter in 
Fig. 4(d). It can be seen that at small Reynolds numbers (roughly 
Reoh < 5000), all data represented tend to collapse on a hydraulic- 
diameter-based correlation even though other shape parameters 
are different. Yet, this tendency is weakened as Re oh becomes 
larger. A least square fit of all data points representing Reynolds 
numbers ranging from 1 to 66,820.6 yields the following correla¬ 
tion, which is shown by the curve in Fig. 4(d). 

Nu Dh = 8.456 + 0.325Re^ 25 Pn (19) 


60000 
50000 
40000 
8-30000 
20000 
10000 
0 

0 10000 20000 

P*Dh 

Fig. 6. Streamwise dispersion conductivity for x-direction flow. 

model for the interrupted plate insert and / = 7.5 mm, t = 0.55 mm 
and 2b = 2.75 mm. The term is: 



In the aforementioned simulation runs, RANS models are used 
for high Reynolds number simulations; the small eddies are rep¬ 
resented by Reynolds stress terms based on closure models. To 
show the pore-scale mixing effect, a Large Eddy simulation (LES) is 
conducted. The Dynamic Smagorinsky-Lilly Model based on [18] is 
used for sub-grid scale modeling. 8,264,256 grid cells are used with 
gradually shrinking cells towards the wall; the maximum y + value 
of the cells adjacent to the wall is 2.59. The time step size is 10 5 s. 
Though the average motion of the fluid is along the plates in the x 
direction, the pore-scale fluid exhibits strong unsteadiness and 
mixing effects, characterized by the cascade of eddies generated as 
a result of the fluid being periodically interrupted by the plates 
(Fig. 5(a)). The maximum local wall heat flux is found at the edges 
of the frontal areas due to impingement of fluid onto these surfaces 
(Fig. 5(b)). The flow is effectively cooled by the interrupted-plate 
medium. 

3.2.3. Thermal dispersion analyses for a fixed shape 

When volume averaging is applied to the energy equation, a 
thermal dispersion term appears. It is associated with spatial var¬ 
iations of velocity and temperature within the pores of the matrix. 
It is anisotropic in that a value found for streamwise dispersion is 
different than a value found for cross-stream dispersion. In this 
section, the streamwise thermal dispersion term for x-direction 
flow (see Fig. 3(c)) is calculated based on CFD runs on the REV 


-pc p {ut) f =-pc PV 2— I (u-<u> f )(r f -<7}) f )dv 

1/rev, f 

( 20 ) 

The streamwise thermal dispersion term is modeled using the 
dispersion conductivity, /cdisp, 

- pCp {ut) f = k disp ^f ( 21 ) 

The dispersion conductivity is normalized based on the air 
conductivity and plotted against the Peclet number. The results are 
shown in Fig. 6. The Peclet number dependency is expected since 
dispersion transport is by mechanisms that are velocity dependent. 
One of the dispersion mechanisms is, by nature, similar to eddy 
transport in turbulent flow. For porous media, the “eddies” are 
within the pores and scale on pore size. 

3.3. Experimental characterization of porous metal foams 

The fine-scale geometry of the metal foam is randomly variable 
(over a feature size range) and essentially unknown to the nu¬ 
merical modeler. Computing the pressure drop, heat transfer and 
dispersion terms for representative, repeated geometries and 



0.003 


0.006 (m) 


Temperature 
fluid 

337 

330 

322 

315 

307 

300 

293 
[K] 



Temperature 
Fluid 

337 
330 
322 
315 
307 

R 

300 

293 
[K] 


o o . 

S' •> 

: lui3’^eamline 
0.003 0.006 (m) 


[m s A -1] 


Wall Heat Flux 
Solid Surface 

-462 
-2632 
-4802 
-6972 
-9142 
-11312 
-13482 
-15652 
-17822 
-19992 
-22162 
-24332 
-26502 
-28672 
-30842 
-33012 
-35182 
-37352 
-39522 
[W m A -2] 


(a) Fluid temperature distribution 


(b) Streamlines, wall heat flux, and fluid temperature distribution 


Fig. 5. Transient flow field of an LES simulation, Re^h = 2872 l = 7.5 mm, t = 0.55 mm, 2b = 2.75 mm, p = 12.05 kg/m 3 . 
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combining them statistically would be approximate and tedious. 
Thus, we have chosen to determine the pressure drop and inter¬ 
facial heat transfer terms through experiments and simulation. The 
dispersion term for the foam matrix is discussed later. 

Pressure resistance characteristics of the metal foams have been 
obtained by measuring pressure drop values across the metal foam 
section when air is blown through the foam at different velocities 
[19]. As shown in Fig. 7, the pressure drop shows a quadratic rela¬ 
tion with Darcian velocity. Data from the measurements lead to the 
coefficients: 

10PPI metal foam: K= 2.397 x 10” 7 m 2 , F= 0.140 

40PPI metal foam: K = 9.913 x 10“ 8 m 2 , F = 0.181 

The indicator PPI refers to the nominal number of pores per inch 
(per 25.4 mm), a frequently-used descriptor of these foams. They 
both have a porosity value of 93%. 

A number of heat transfer correlations for porous media are 
available from the literature. To choose one that is most suitable, 
another experiment is done 19]. A schematic of the experimental 
setup is shown in Fig. 8. Water is circulated in the clockwise di¬ 
rection of the figure and pumped into the lower section of the 
liquid-piston compression chamber to compress air. Two experi¬ 
mental runs each are conducted with 10PPI and 40PPI metal foam 
inserts. Instantaneous values of volume and pressure are measured 
and instantaneous gas-volume-average temperatures are calcu¬ 
lated using the ideal gas law. The experimental runs were simu¬ 
lated using the zero-dimensional model by substituting different 
heat transfer correlations found in the literature and the results 
were compared with the experimental results, as shown in Fig. 9. A 
modified version of the Kamiuto heat transfer correlation was 
selected, as discussed in Ref. [19 . 

Il fm = 0.996 ( P u °* dm p r \ (e = 0.93) (22) 

where Uoa is the Darcian velocity, the velocity that the fluid would 
have in the absence of the solid. 


3.4. Simulations of a compressor filled with metal foam 


Numerical modeling of the metal-foam-inserted liquid-piston 
compressor combines Volume of Fluid (VOF) and two energy 
equations, one for the fluid and a second for the solid. Thus, we 
solve for energy transport in the fluid and solid phases and couple 
the two through the interfacial heat transfer correlation. The VOF 
method tracks the instantaneous location of the water/air interface 
in the compression chamber. Let subscripts 1 and 2 represent air 
and water, respectively. The continuity equations for air and water 
are: 
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Fig. 7. Pressure drop measured from the experiments [19]. 
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Fig. 8. Schematic of the experimental setup used to select the optimum open-cell 
metal foam heat transfer correlation [19 . 


^l^ + V-L( Pl ) f (u) f j =0 (23) 

and 

^+V-(a 2 <u> f ) =0 (24) 

Momentum and energy equations are solved for the immiscible 
fluid mixture. 

a<p) a f ) U>f + V- (<p} f <u) f <u) f ) = — V(p) f + V-(7> f + (p) f g + S m 

(25) 
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Fig. 9. Comparison of experimental results to the Zero-D model solutions obtained by 
using different heat transfer correlations [19]. (Variables are non-dimensionalized by 
initial values). Correlations: Fu et al. [20 ], Kamiuto and Yee 21 , Kuwahara et al. [22], 
Nakayama et al. [23], Wakao and Kaguei 24], and Zukauskas.[25j. 
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where (p) f = a 1 (p 1 ) f + a 2 {P2) f (26) 

The stress tensor, ¥, is based on viscosity of the mixture of air 
and water in nodes that contain the liquid-gas interface. 


\i = aq/ii + a 2 M2 (27) 

The momentum sink term, S m . x = 9p/9x, is given by Eq. (12). 
The energy equation for the fluid is. 


0 £ (c p <p) f <r> r 

dz 


+ eV-({u) f (p) f C p (T) r 


= eW-kfW{T) f — V-pc p { u f> f + h v ((T s ) s 



+ £ 


9(P) f 

6/ 


(28) 


where 


(p) f C p = a\ (pi) f Cp 7l + ^2(P2) fc p,2 

(29) 

and 


kf = + 0 ( 2^2 

(30) 


The energy equation for the solid is. 


(1 - e)^((p s ) 5 Cs(T s ) s ) = (1 - £)V-/< s V(Ts} s - h v ((T s > s - <7} f ) 

(31) 

In the study of ref. [19 , thermal dispersion term was neglected 
(Vpc p (u f ) f = 0). It shows that when only a portion of the chamber 
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Fig. 10. CFD simulations on the chamber with metal foam insert (note that up is to the right in this figure). 10PPI metal foam. U 0 = 0.206 m/s, tf= 1.3 s. The chamber length and 
radius are respectively: 0.294 m and 0.0254 m. (Shown on the left is the entire chamber region; shown on the right is a scaled-up view of the air region). 
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is occupied by porous insert, flow is stable in the porous insert 
region but vortices develop elsewhere (Fig. 10(a)); when the entire 
length of the chamber is occupied by insert, velocity streamlines 
are smooth (see Fig. 10(b)). This indicates that optimization of the 
distribution of porous inserts in the chamber may be done to 
achieve minimal temperature rise of air during compression, and 
best compression efficiency. 

In the present study, a CFD simulation is conducted with the 
thermal dispersion term included. The chamber is fully occupied by 
10PPI metal foam. The thermal dispersion term is modeled based 
on dispersion conductivity. 

-pc p <nT) f = K dis V<7) (32) 

Kuwahara and Nakayama [26] proposed a closure model for 
thermal dispersion for 2-D cross flow over small square rods: 


kdis,XX 



kdis,YY 



kdis,XX 



kdis,YY 



p e 2 

0 022 <,- D e) K<’°) 

(33) 

Pe 17 / \ 

0.022 d (Pe D < 10) 

(1 — e) 1 ' 4 V > 

(34) 

K>io) 

(35) 

0.052(1 - <?) 1/2 Pe D (Pe D > lo) 

(36) 


where XX is the local streamwise direction and YY is the local cross¬ 
stream direction. This model is implemented in the CFD simulation 
in the present study. The result shows that adding thermal 
dispersion to the energy equation does not affect the flow field 
much, as shown in Fig. 10(c). Close scrutiny shows that at the end of 
compression, the case with dispersion included has a smaller 
maximum local temperature than the case neglecting dispersion 
Transient values of bulk temperature of air are compared be¬ 
tween simulations with and without the dispersion term (Fig. 11). 
Results of both cases show that, though the compression time is 
short and isentropic compression would lead to a temperature of 
575 K, the insert is effective in suppressing temperature rise. 
Further, they show that during the latter period of compression, 
calculated temperatures of the case with dispersion are less, 
because of larger thermal gradients and higher heat transfer. 


Thermal dispersion has a greater influence on the maximum bulk 
temperature when temperatures and temperature gradients rise 
near the end of compression, particularly for large pressure ratios. 
When thermal dispersion effect is included, the computed bulk air 
temperature at the final compression state is 5 K less than when the 
dispersion effect is not included. The local maximum temperature 
values at the final compression state for the case without dispersion 
effect and the one with dispersion effect are respectively: 417.2 K 
and 428.3 K. This temperature difference caused by dispersion has a 
small effect on the compressor performance for the cases investi¬ 
gated in the present study, since this temperature difference is 
small compared to the total amount of temperature rise during 
compression. 

4. Conclusions 

A compressed air energy storage system designed for use in 
wind turbine plants was introduced and the importance of thermal 
control during compression and expansion was emphasized. A 
compressor design was proposed and models for its analysis were 
presented. Noted in modeling was a need for sub-models to char¬ 
acterize pressure drop in the porous insert, heat transfer from fluid 
to solid phases of the insert and thermal dispersion within the 
insert. Such models were presented for two candidate inserts, one 
based on interrupted plates and another based upon metal foams. 
The models came from detailed computation of a representative 
elementary volume, in the case of the interrupted plate, and from 
experiments, in the case of the foam. Computational results using 
these models for the foam insert case were presented. The inserts 
are very effective in suppressing temperature rise. The metal foam 
kept the bulk temperature below 360 K if it occupied the full 
chamber length; and, for reference, the adiabatic bulk temperature 
would have been 575 K if no heat transfer medium had been 
introduced. Computed results that compare a case with the insert 
occupying the full chamber length and a case having the insert 
occupying only a portion of the chamber length show that the 
insert stabilizes the flow, where present; but beyond the insert 
there is a tendency for development of secondary flows. The 
diameter of the secondary flow vortex scales on the thickness of the 
non-porous region. A study in which dispersion was added to the 
analysis showed that effect of dispersion is small and in the di¬ 
rection of reducing the maximum average temperature reached 
during compression. Dispersion weakens temperature gradients in 
the fluid, increasing near-wall gradients and heat transfer to the 
walls. For the pressure ratio studied, dispersion has reduced the 
maximum local temperature by 10 K. 


360.0 



Time (s) 

Fig. 11. Bulk temperature of air. 
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